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ABSTRACT 

We present a measurement of the evolution of the stellar mass function (MF) of galaxies and the evolution 
of the total stellar mass density at < z < 5, extending previous measurements to higher redshift and fainter 
magnitudes (and lower masses). We use deep multicolor data in the Fors Deep Field (FDF; /-selected reaching 
Iab ~ 26.8) and the GOODS-S/CDFS region (/T-selected reaching Kj^^ ^ 25.4) to estimate stellar masses based 
on fits to composite stellar population models for 5557 and 3367 sources, respectively. The MF of objects 
from the /T-selected GOODS-S sample is very similar to that of the /-selected FDF down to the completeness 
limit of the GOODS-S sample. Near-IR selected surveys hence detect the more massive objects of the same 
principal population as do /-selected surveys. We find that the most massive galaxies harbor the oldest stellar 
populations at all redshifts. At low z, our MF follows the local MF very well, extending the local MF down 
by a decade to 10** M0. Furthermore, the faint end slope is consistent with the local value of a ^ 1.1 at least 
up to z ^ 1.5. Our MF also agrees very well with the MUNICS and K20 results at z < 2. The MF seems to 
evolve in a regular way at least up to z 2 with the normalization decreasing by 50% to z = 1 and by 70% to 
z = 2. Objects with M > 10"^ M© which are the likely progenitors of todays L> L* galaxies are found in much 
smaller numbers above z ^ 2. However, we note that massive galaxies with M > 10" M0are present even to 
the largest redshift we probe. Beyond z ~ 2 the evolution of the mass function becomes more rapid. We find 
that the total stellar mass density at z = 1 is 50% of the local value. At z = 2, 25% of the local mass density is 
assembled, and at z = 3 and z = 5 we find that at least 15% and 5% of the mass in stars is in place, respectively. 
The number density of galaxies with M > lO" Mq evolves very similarly to the evolution at lower masses. It 
decreases by 0.4 dex to z ^ 1, by 0.6 dex to z ~ 2, and by 1 dex to z ^ 4. 

Subject headings: surveys — cosmology: observations — galaxies: mass function — galaxies: evolution — 
galaxies: fundamental parameters 



1. INTRODUCTION 

The stellar mass of galaxies at the present epoch and the 
build-up of stellar mass over cosmic time has become the fo- 
cus of intense research in the past few years. 

Generally, this kind of work relies on fits of multi-color 
photometry to a grid of composite stellar population (CSP) 
models to determine a stellar mass-to-light ratio, since large 
and complete spectroscopic samples of galaxies (at z > 0) are 
not yet available. 

In the local universe, results on the stellar mass func- 
tion (MF) of galaxies were published using the new gen- 
eration of wide-angle surveys in the optical (Sloan Digital 
Sky Survey; SDSS. York et al. 2000: 2dF, Colless et al. 2001) 
and near-infrared (Two Micron All Sk y Survey; 2MASS, 
iSkrutskie et alJl997h . ICoieetal] ll2()0l combined data from 
2MASS and 2dF to derive the local stellar MF, Bell et afl 
J2003h used the SDSS and 2MASS to the same end. 

At z > 0, a number of authors studied the stellar 
mass density as a function of redshift ( Brinchmann & Ellip 
|2p00* 'Droryetal."'2001', "Cohen"'2002» 'Dicki nson et aL.2003r 
iFontanaet al. 2003; Rudnicket al. 2003) reaching z ~ 3. It 
appears that by z ^ 3, about 30% of the local stellar mass den- 
sity has been assembled in galaxies, and at z ^ 1, roughly half 
of the local stellar mass density is in place. This seems to be 
in broad agreement with measurements of the star formation 
rate density over the same redshift range. 
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Others investigated the evolution of the MF of galaxies 
(iDror v et alJl2004t iFontana et alJl2004^ to z 1.5, finding a 
similar decline in the normaUzation of the MF. However, it 
is possible that galaxies evolve differently in number density 
depending on their morphology. 

In this letter, we extend the measurement of the stellar MF 
and its integral, the total stellar mass density, to z ~ 5. We 
describe the data set in Sect. Hand present the method used to 
derive stellar masses in Sect. |3l We discuss the results on the 
stellar MF in Sect.|3 Finally, we discuss the total stellar mass 
density and the number density of massive galaxies in Sect.|5l 

Throughout this work we assume ^Im = 0.3, ft\ = 0.7, //o = 
70 km s~' Mpc"' . Magnitudes are given in the AB system. 

2. THE GALAXY SAMPLE 

This work is based on data from two deep field surveys hav- 
ing multicolor photometry and rich foUo wup spectroscopy, 
the Fors Deep Field (FDF; Heidt et al.l l 2003l iNoU et alJ 
|2004) covering the UBgRIzJKs bands in 40 arcmin and the 
GOODS-S/CDFS field (|(Jiavalisco et al. 2004') covering the 
UBVRIJHKs b ands in 50 arc min-. This sample is identical to 
the one used in lGabasch et al.l (l2004bl) to study the global star 
formation rate to z 5. 

The FDF photometric catalog is published in iHeidt et all 
("2003") and we use the I-band selected subsample cover- 
ing the deepest central region of the field as described in 
iGabasch et al.l ll2004a). This catalog lists 5557 galaxies down 
to / ~ 26.8. The latter also discusses the I-band selec- 
tion and shows that this catalog misses at most 10% of the 
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objects found in ultra-deep K-band observations found by 
iLabbe e t al. (2003). 

Photometric redshifts for the FDF are calibrated against 
362 spectroscopic redshifts up to z ~ 5 and have a accu- 
racy of lS.z/{Zspec + 1) < 0.03 with only ~ 1% outliers. The 
method and this calib ration are presented and discussed in 
lGabaschetan(l2004at . 

Our K-band selected catalog for the GOODS-S/CDFS field 
is based on the publicly available 8 2.5 x 2.5 arcmin^ J, H, Ks 
ISAAC images, taken at ESOA'LT with seeing in the range 
0.4" -0.5". The U and I images are from ESO GOODS/EIS 
public survey, while the B V R images are taken from the 
Garching-Bonn Deep Surv ey. Th ese datasets are e xtensivel 
discussed in Arnouts et al] j2001h and lSchirmer et al. (2003| 
respectively. 

The data for the GOODS field were analyzed in a very sim- 
ilar way to the data of the FDF. The objects were detected in 
the K-band images closely following the procedure used for 
the FDF I-band detection. A detailed description of the pro- 
cedure as well as the catalogs can be found in Salvato et al. 
(2004, in preparation). We detect 3367 sources in the K-band 
for which we measure fixed aperture and Kjon-like total mag- 
nitudes in all bands. Of these, 2 objects are detected only in 
K, and 36 objects are detected only in the Ks, H, and J band 
images. Number counts match the literature values down to 
K K, 25.4, which is the completeness limit of the catalog. We 
computed photometric redshifts using the same method and 
SED template spectra as for the FDF. A comparison with the 
spectroscopic redshifts of the VIMOS team and the F0RS2 
spectra shows that the photometric redshifts have an accuracy 
Az/fcp,,+ l)< 0.056. 

3. DERIVING STELLAR MASSES 

We derive stellar masses (the mass locked up in stars) by 
comparing multi-color photometry to a grid of stellar popula- 
tion synthesis models covering a wide range in parameters, 
especially star formation histories (SFHs). This method is 
described in detail and compared against spectroscopic and 
dynamical mass estimates in Drory, Bender, & Hopp (2004), 
hence we only briefly outline the method her e. 

We base our model grid on the Bruzual & Charlol ll2003h 
models. We parameterize the possible SFHs by a two- 
component model, a main component with a smooth ana- 
lytical SFH modulated by a burst of star formation. The 
main component is parameterized by a SFH of the form 
ipit) oc exp(-f/T), with r G [0.1, cxd] Gyr and a metallicity of 
-0.6 < [Fe/H] < 0.3. The age, f, (defined as the time since 
star formation began) is allowed to vary between 0.5 Gyr and 
the age of the universe at the object's redshift. This compo- 
nent is linearly combined with a burst modeled as a 100 Myr 
old constant star formation rate episode of solar metallicity. 
We restrict the burst fraction, (3, to the range < /3 < 0. 15 in 
mass (higher values of /3 are degenerate and unnecessary since 
this case is covered by models with a young main compo- 
nent). We adopt a Salpeter IMF truncated at 0.1 and 100 
for both components. A different slope will change only the 
mass scale. However, if the IMF depends on the mode of star 
formation the shape of the MF can be affected. 

Additionally, both components are allowed to exhibit a dif- 
ferent and variable amount of extinction by dust. This takes 
into account the fact that young stars are found in dusty envi- 
ronments and that the starlight from the galaxy as a whole may 
be r eddened by a (geometr y dependent) different amount. In 
fact. lStasihska et all (I2004^ find that the B aimer decrement in 
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Fig. 1 . — Stellai' masses vs. redshift in tlie FDF (squaies) and in GOODS- 
S (triangles). Tlie luminosity-weiglited mean ages of the stellar populations 
as determined by the CSP fit are coded by colors ranging from red for old 
populations (age > 5 Gyr) to blue for young populations (age < 1.5 Gyr). 



the SDSS sample is independent of inclination, w hich, on the 
other hand, is driving global extinction (see, e.g.. lTullv et alJ 
1998). 

We compute the full likelihood distribution on a grid in 
this 6-dimensional parameter space (r, [Fe/H],f,Ay,/3,Ay), 
the likelihood of each model being cx exp(-x^/2). To com- 
pute the likelihood distribution ofM/L, we weight the M/L of 
each model by its likelihood and marginalize over all parame- 
ters. The uncertainty in M/L is obtained from the width of this 
distribution, which is on average between ±0.1 and ±0.3 dex 
at 68% confidence level. The uncertainty in mass has a weak 
dependence on mass (increasing with lower S/N photome- 
try) and much of the variation is in spectral type: early-type 
galaxies have more tightly constrained masses than late types. 
At high redshift, however, the uncertainty in M/L grows as 
objects drop-out in the blue bands and stellar populations be- 
come younger: the typical uncertainty at z 3 is 0.6 dex and 
becomes as large as 1 dex at the highest redshifts we probe. 

In Fig^we show the distribution of galaxies in the mass 
vs. redshift plane for the FDF (squares) and GOODS-S (trian- 
gles). In addition, we code the age of each galaxy (using the 
best-fitting model) in colors ranging from blue for young (age 
< 1 .5 Gyr) to red for old stellar populations (age > 5 Gyr). 

The distribution of objects from the /T-selected GOODS- 
S sample is very similar to the distribution of the /-selected 
FDF down to the completeness limit of the GOODS-S sample, 
which is ^ 1 dex shallower in limiting mass. This indicates 
that present optical and near-IR surveys are unlikely to have 
missed a substantial population of massive objects, with the 
possible exception of heavily dust-enshrouded sources which 
may escape detection in both optical and near-IR surveys. 

A striking feature of Figgis that the most massive galax- 
ies harbor the oldest stellar populations at all redshifts. There 
always exist galaxies which are old (relative to the age of the 
universe at each redshift) but less massive, yet the most mas- 
sive objects at each redshift are never among the youngest. 



Stellar mass function at < z < 5 
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Fig. 2. — The stellar mass function (MF) as a function of redshift. The local MF of Bell et al. 1 2003) and Cole et al.' f200H are shown in all redshift bins for 
comparison. Results from the MUNICS survey tProrv et al-2004.) and the K20 survey iFontana et al-2004J at z < 2 are also shown. 



4. THE STELLAR MASS FUNCTION 

Fig.|2]shows the Vmax-weighted mass function in seven red- 
shift bins from 2 = 0.25 to z = 5. For comparison, we also show 
the local mass function ( Bell e t al- 2 00 3: Cole et al. 20Ql) and 
the MPs to z ~ 1.2 of MUNICS jProrv et alJl20dir and to 
z 2 by the K20 survey ( Fontana et alJ2004 . 

The MF of objects from the /T-selected GOODS-S sample 
is very similar to that of the /-selected FDF down to the com- 
pleteness limit of the GOODS-S sample (with the exception 
of a possible slightly lower normalization at z > 2.5 by about 
10% to 20%). This shows that near-IR selected surveys at 
high redshift essentially detect the more massive objects of 
the same principal population as do optically (/-band) selected 
surveys (see also Sect. 3 and Fig. 1 in Gabasch et al. 2004aV 
It remains to be seen what fraction of massive galaxies (and 
total stellar mass density) might be hidden in dusty starbursts 
which appear as sub-mm galaxies. 

In our lowest redshift bin, z ^ 0.5, the MF follows the local 
MF very well. The depth of the FDF (/ - 26.8) allows us to 
extend the faint end o f the MF down t o 10^ M©, a decade 
lower in mass than in iBell etalJll2?OTl . with no change of 
slope. Furthermore, the faint end slope is consistent with the 
local value of a 1.1 at least to z ^ 1.5. Our mass function 
also agrees very well with the MUNICS and K20 results at 

The MF seems to evolve in a regular way at least up to 
Z ~ 2 with the normalization decreasing by 50% to z = 1 and 
by 70% to z = 2, with the largest change occurring at masses of 
M 10"^ M0. These likely progenitors of todays L> L* galax- 
ies are found in much smaller numbers above z ~ 2. However, 
we note that massive galaxies with M > 10" M0are present 
even to the largest redshift we probe (albeit in smaller num- 
bers). Beyond z 2 the evolution becomes more rapid. 



It is hard to say whether the difference between the FDF and 
the GOODS-S field visible a t z 3 in Fig.|2lis significant. In 
fact, iGabaschetall ilOOm find not much difference in the 
rest-frame UV luminosity function in the very same dataset. 
We think it might be related to the shallower near-IR data in 
the FDF compared to GOODS-S. Less reliable information 
on the rest-frame optical colors at young mean ages might in 
fact lead to an overestimation of the stellar mass. This would 
also explain why the two MFs become similar again at even 
higher z when the near-IR data in GOODS-S reach their limit, 
too. This would mean that the mass density in the FDF at 
z > 2.5 might be overestimated (see below). 

5. STELLAR MASS DENSITY AND NUMBER DENSITIES 

Fig. |3] shows the evolution of the total stellar mass den- 
sity along with data from the literature, extending the avaiL 
able data to z ~ 5 ( Coleetal. 2001; Brinchmann & Elli 

2000; Drorvetal.ll2()'0 1: Cohe a.2002: .Dickinson et al 2003*1 

Fonta na et al.l2003L 12004; Dro rv et alJ2004 . 

We compute the total stellar mass density by directly sum- 
ming up contributions from all objects in both fields (we ob- 
tain very similar values by means of fitting Schechter func- 
tions to the data in Fig. |2}. We find that the stellar mass 
density at z = 1 i s 50% of the local value as determined by 
'Cole et alJ j2001l) . At z = 2, 25% of the local mass density is 
assembled, and at z = 3 and z = 5 we find that at least 15% and 
5% of the mass in stars is in place, respectively. Fig.|3lshows 
the results from both fields separately, however 

We als o show the in t egral o f the star formation rate deter- 
mined bv lSteidel et al.l (11999') as a dotted line. Furthermore, 
the dashed line shows the same quantity determined from the 
rest-frame UV luminosity function of in same dataset used 
here ( Gabasch et al. 2004b ). We find these measurements in 
agreement with each other, and with the mass densities de- 
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Fig. 3. — The total stellar mass density of the universe as a function of 
redshift. We also show the integral of the star formation rate as a dotted line 
and dashed line. 
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Fig. 4. — The number density of galaxies with stellar masses M > IO'^Mq 
and M > 1O"M0. Results from the MUNICS survey as well as results from 
Sub-mm studies are also shown. 



rived here and in the literature before. However, the mass 
densities do show considerable scatter especially at redshifts 
above 1 .5. However, as discussed above, the mass density 
in the FDF might be overestimated, hence reducing the scatter 
between our two fields. The FDF also lies above the average 
literature values in the redshift range in question, while the 
GOODS-S data seem to be in better agreement with previous 
measurements. 

The fact that the stellar mass densities and the integrated 
star formation rates give consistent results is encouraging. 
However, dust-enshrouded starbursts at high redshift would 
be missing from both the SFR and the MF. However, if 
the number density of massive sub-mm galaxies derived by 
iTeczaet al.l (l 2()()l of 8.9!^J x IQ-^Mpc"^ at 1.8 < z < 3.6 is 
typical of these systems, they do not contribute significantly 
to the mass density at these redshifts (see below and Figl^i. 

Finally, Fig.|4]addresses the number density of massive sys- 
tems as a function of redshift. We show the numbers of sys- 
tems with M > 1O'"M0 and M > 10" Mq as full and open 
symbols, respectively. The number density of massive sub- 
mm galaxies estimated by Tecza et al. (.2004 ) is also shown. 
The most striking features of Fig. 21 is that the number den- 
sity of the most massive systems shows evolution which is 
very similar to the evolution of the number density at lower 
masses over this very large redhsift range. Massive systems 
are present at all redshifts we probe, their number density de- 
creasing by 0.4 ±0.2 dex to z 1, by 0.6 ±0.3 dex to z 2, 
and by 1 ± 0.45 dex to z 4. 
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